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Abstract 

A Swinging Atwood Machine {SAM) is built and some experimental results concerning 
its dynamic behaviour are presented. Experiments clearly show that pulleys play a role 
in the motion of the pendulum, since they can rotate and have non-negligible radii and 
masses. Equations of motion must therefore take into account the inertial momentum 
of the pulleys, as well as the winding of the rope around them. Their influence is com- 
pared to previous studies. A preliminary discussion of the role of dissipation is included. 
The theoretical behaviour of the system with pulleys is illustrated numerically, and the 
relevance of different parameters is highlighted. Finally, the integrability of the dynamic 
system is studied, the main result being that the Machine with pulleys is non-integrable. 
The status of the results on integrability of the pulley-less Machine is also recalled. 

Key words: Swinging Atwood's Machine {SAM), Chaotic system. Nonlinear dynamics, 
Integrability, Experimental SAM apparatus. 
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1. Introduction 

This paper deals with the Swinging Atwood Machine [SAM), a non-hnear two- 
degrees-of-freedom system derived from the well-known simple Atwood machine. The 
latter was devised in 1784 by George Atwood, a London Physics lecturer who built his 
own apparata as a means of practical illustration, in order to experimentally demonstrate 
the uniformly accelerated motion of a system falling under the earth gravity field g with 
mass dependence In Atwood's original machine, two masses are mechanically linked 
by an inextensible thread wound round a pulley. In SAM^ one of the masses (m) is 
allowed to swing in a plane while the other mass (M) plays the role of a counterweight; 
hence SAM can be seen as a parametric pendulum whose length is varying as a function 
of the parameter /i = M/m. 

For about twenty-five years, many studies have been conducted concerning the me- 
chanical behaviour of SAM. Said studies were conducted exclusively on a simplified model 
for SAM neglecting any influence from a massive set of pulleys. Through numerical inves- 
tigations, [30| inferred the pulley-less SAM to be an extremely intricate system exhibiting 
significant changes in the qualitative behaviour of trajectories, depending on /z-values. 
Assuming > 1, motion is limited in space and two types of trajectories can be dis- 
tinguished owing to the initial conditions: singular ones for which pendulum length is 
initially zero, and non-singular ones where the pendulum is initially released from rest 
with a non-zero length. For the former, it appears that /i = 3 is a particular condition 
corresponding to terminating trajectories, i.e. those for which pen dulum length becomes 
zero after a given duration, regardless of the initial conditions [29| . The latter is divided 
into periodic, quasi-periodic and what can be conjectured to be ergodic trajectories in 
some domain. SAM without massive pulleys was also studied by means of Poincare sec- 
tions wherein chaotic dynamic behaviour becomes prominent as /i is increased [27| . An 
interesting and surprising result is the integrability of the pulley-less SAM for /Lt = 3, a 
conclusion which is also supported using Hamilton- Jacobi theory [2^ and Noether sym- 
metries [19|. For /X > 3, [§| proved that SAM without massive pulleys is not integrable, 
contrary to what was speculated by [23|. The belonging of /i = M/m to a special set of 
parameters {/ip : p G Z} was established as a necessary condition for integrability; this 
result was proven independently in [§], 35 1 and 0, and is proven in Remark 17. If 2) of 



the present paper as well. Moreover, unbounded trajectories < 1) have been studied 
via energetic considerations [sij; [l^l identified and classified all periodic trajectories in 
the pulley-less SAM for ^ = 3. Finally, a very recent result co-written by one of the 
authors of the present article ([31) proved the non-integrability of this pulley- less model 
for SAM for the exceptional values {/ip : p G Z}; this had been an open problem, at least, 
since [Ol] explicitly tackled the issue for the first time. It is worth noting that all of these 
studies are theoretical, albeit for the most part strongly supported by massive numerical 
simulations. 

In this paper, we intend to describe a useful physical construction of SAM in detail, 
as well as present further experimental and theoretical results. In addition, a theoretical 
premise is introduced which stands as a novelty all its own: to wit, as suggested by 
experiments, pulleys are no longer neglected, in order to take account of non-zero radii 
and rotation around their axes of revolution. When dealing with 7V-degree-of-freedom 
non-linear systems, the modern researcher's tendency to restrict adjectives such as "com- 
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plex" to ^ 1 should not divert us from the fact that even the dynamics of apparently 
elementary cases such as = 2 are often very difficult to determine {^), and thus a 
source of interest in their own right. As shown in this paper, such is the case for SAM. 
A schematic representation of SAM is featured in Section [5] partly aimed at the deriva- 
tion of the equations of motion with pulleys in Section [31 The constructed apparatus 
is then described in detail in Section [H and some experimental results are presented in 
Section [5l A comparison with the theoretical model is performed in Section [6] through 
numerical simulations of the general equation of motion obtained in Section [31 Section [71 
is devoted to prove the non-integrability for SAM. The rigorous proof shown therein is 
requisite to establish definitively the non-integrability of SAM, hence incumbent upon 
any proper completion of experimental and numerical results. Indeed, although for some 
systems non-integrability is somehow suggested by a thorough Poincare section anal- 
ysis, chaotic zone detection may require extremely careful numerics and will at times 
become laborious - we might call this shy chaos. Furthermore the lack of integrability 
of some systems cannot be discovered by looking at the real phase space. There exist 
non-integrable systems without any recurrent motion in the real phase space and such 
that chaos is confined away from the real domain. For an example, see [l8l |. Finally, 
Section [8l concludes with perspectives on further experiments and some comments on 
results concerning integrability of the pulley-less case. 



[SO'I suggested a SAM physical demonstration model using a vertically mounted air 
table, and alleged a successful experimental demonstration of the system's motions. How- 
ever, there is no experimental result in the aforementioned reference and its proposed 
model for SAM is unequivocally not as close to the theoretical system as the model 
described herein (cf. Section[21). To our present knowledge, detailed experimental studies 
of SAM, let alone comparisons of any such experiments with the theory, do not exist in 
the literature prior to our work. 

Therefore, our work arguably completes the above theoretical and experimental re- 
search on SAM, and, at the same time, opens new problems and sets a starting point for 
further experimental and theoretical studies. 

2. Schematic representation of SAM 

SAM is represented by the system S sketched in Figure [Hand consisting of: 

• a pendulum, considered as a material point A of mass m, 

• a counterweight, considered as a material point B of mass M, 

• a thread of length L linking A and B, 

• two pulleys V and V' of radius R, distant from one another by a predetermined 
distance D. 

S is studied relative to the Galilean laboratory frame TZ — {O Cy e^) whose origin 
O is chosen to correspond to the centre of the pendulum pulley. Axis Ox corresponds 
to the pulley revolution axis; axis Oy is the horizontal direction defined by the pulley 
centres (O and O'), and oriented from O toward O'; finally, Oz is chosen to correspond. 
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Figure 1: Schematisation of the system S = {pendulum-thread-pulleys-counterweight} representing 
SAM. The angle (p locates the material point N' of the pulley T" ■ 



for the sake of convenience, to the downward direction of the local earth gravity field g 
(vertical). 

Pendulum A is characterized by its variable length r = QA, Q being the geometrical 
point where the thread departs from the pulley, and by the angle 9 formed by QA and 
the downward vertical. Note that 6 as represented in Figure [1] is a positive angle. 

Vertical motion of the counterweight B is described by its coordinate zb, which can 
be related to the angular position ip of any point N on the pulleys, provided the thread 
does not slip on the pulley (in Figure (TJ for the sake of clarity, N is drawn on the pulley 
V' associated to B and thus labelled N'). Indeed, under this assumption, when B is 
falling down, pulleys are able to rotate in such a way that the velocity of any point of 
the pulleys (for instance N') is equal to the velocity of B. Hence: 

zb = R(p (1) 

Note the difference between the rotation angle ip of the pulleys and 9: the former defines 
the location of any material point on a pulley, whereas the latter defines the angular 
position of A, as well as that of the geometrical point of contact Q. This subtlety is due 
to the necessary mechanical description of contact in terms of three points [2]| : 

• the geometrical point of contact Q, 

• the point N of the pulley V and corresponding to Q at time t, 

• the point K of the thread corresponding to Q at time t. 

Figure m displays the configuration. A similar problem can be found in [2^ . 



Figure 2: Description of the contact between a pulley and the thread with three points: K, Q, and N. 
At any time t, both K and TV touch each other at Q. In the figure, these points are separated for the 
sake of clarity. 

A physical way of understanding the difference between 6 and ip is to imagine the 
following situation. At initial time, assume that Q and A'^ are superposed: 9 = 9o and 
ip = (fio (Figure [3^). If 9 is fixed and B heads downwards with velocity zb, the absence 
of slippage of the thread on the pulleys implies that they rotate with angular velocity tp 
given by ([T]), meaning A" is moving while Q remains fixed, and r evolves from rg to r; at 
final time, ip tpo (Figure Os). 



Figure 3: Evolution of point A'^ of the pulley V when is kept constant. The geometrical point of contact 
Q is immobile and only r evolves. 

Since L is supposed to be constant, zb is directly related to r, D, and the lengths 
QP and P'Q' corresponding to regions where the thread and the pulleys keep contact. 
Because 
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one has, precisely, 

L = r + D + iiR- RO + ZB- 
Since — r — R6 + zb, it follows that 

. _ ZB _ i?0 - r 
^ "R ^ R 

Hence, 5 is a system with two degrees of freedom, for instance 6 and r. 



(2) 
(3) 



3. Equations of motion of SAM with pulleys 

3.1. Equations of motion 

Let us determine the equations of motion for SAM by taking into account the pulleys, 
as opposed to what has been assumed in previous theoretical studies 1^, 2^ 27, H, 2^, 
Isol. [sil. [35II . Indeed, their non-zero radii imply a likely change in position for Q, and the 
equally likely rotation of V and V' around their respective revolution axes apparently 
deems their inertial momentum Ip a significant dynamic parameter. Observations will 
confirm this - see Section [5j 

Lagrange's formalism is used to derive these equations. The kinetic energy of the 
system is expressed by: 
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The first term sums up the contribution by pendulum A, the second is relative to the 
counterweight B and the third one corresponds to the rotation of the two pulleys. We 
have: 

dOA dOQ dQA 

where OQ — —R cos 9 By — Rs'mOez and QA — — rsinf e 
Cartesian base {O ey e^) we can write 



y 



rcosOe-,, hence in the 



( ° ^ 

va — RS sin 6 — rO cos 9 — r sin 9 

\ -R9cos9 -r9sm9 + rcos9 J 
Similarly, Vb — dOB/dt with OB ~ zb^z- Using (O, one gets 

vb ^ ZBez = - e-z- 
Finally, using ^ and introducing the effective total mass of the system 

2h 



we get 



Mt = M + ni 



£k = ^Mt{R9 - rf 



i?2 ' 



—mr 
2 



This expression is similar to that obtained when neglecting the pulleys, except that: 
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• the total mass is now different from M + m by the term 2Ip/ conveying the 

rotation of the pulleys; 

• the counterweight influences pendulum A through its length r and the winding of 
the rope on the associated pulley. The latter influence is considered in the term 

m. 

Potential energy is only due to the Earth's local gravity field. Dropping an irrelevant 
additional constant term, we have: 

£p,g = —mg ■ OA — Mg ■ OB = —mgzA — Mgzs, 

so that 

£p g = mg{Rsm9 - r cos9) + Mg{r - R9). 
The Lagrangian L(r, 6, f, 6) = £k — £p,g of the system is thus: 

L{r,e,r,9) = ^Mt{R9 ~ r)^ + ^mr^ 9^ - gr{M - m cos 9) - gR{msm9 - M9), 

from which one deduces the conjugate momenta Pr and pe associated to r and 9 respec- 
tively: 

dh 

Pr = -gT = -Mt{R9 - f) 
dL 

pe = ^ = MtR{R9 -f) + mr^9 = -Rpr + mr^9 

89 

The SAM Hamiltonian is, in this case, by definition: 
or else, expressed in terms of Pr and p$: 



H{r,9,pr,pe) = ^ 



p2 {pe + Rpr)"^ 



+ gr{M -mcos9)+gR{msm9~ M9). (4) 



Equations of motion follow from the Hamilton's equations: 

dH , . dH 

Pr = -^ and Pe = -^, 

yielding 

Ht{f-R9) = r9^ + g {cos 9 - n) 



r9 = -2f9 + R9'^ -gsm.9 



with fi = M/m and fit = Mt/m = 1 + /x + (2Ip/mR^). 
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3.2. A more physical way to obtain equations of motion 

An alternative method can be used in order to derive equations of motion ([5]) . The 
second of these is obtained by applying the angular momentum theorem at the mobile 
point Q in order to cancel reaction force at this contact point [2l| : 

— — +vq X mvA = QA X mg, 
at 

where Lq = QA x tjiva = mr^d and 

-Q = JqAQ, -a = ^oq + —qa. 

We obtain r9 — ~2r9 + R9^ — gsm9. Taking this into account, the first equation in ([5]) 
comes from the conservation of mechanical energy in SAM: 

= ^Mt{R9 - ff + ^mr^9^ + mg {R sm9 - r cos 9) ~ Mg{R9 ~ r) = C 
C being a real constant, after derivation with respect to t. 

3.3. Comparison with previous studies 

Without any pulley influence, i.e. pulley inertial momentum /p = and pulley radius 
i? = 0, we recover the equations obtained by [30l |: 

(l + ^)f = r9'^ + g{cos9 - fi) , , 

r9 = -2i-9-gsm9 ^ ^ 

Obviously, if there is no oscillatory motion (9 = 0), the well-known simple Atwood 
machine [5|] is recovered: 

Mtf = g{m - M) (7) 



4. Description of the experimental apparatus 

A physical prototype for SAM has been built using two identical pulleys, a nylon 
thread, a brass ball as a pendulum and a set of different hook masses acting as counter- 
weights. The pendulum and the chosen counterweight are linked together by the nylon 
thread placed around the pulleys. A photo of SAM is displayed in Figured 

4.1. About the pendulum and the counterweight 

Each mass in the experimental device has been measured with a precision scale of 
0.01 g of accuracy. The pendulum is a brass ball with a 30 mm diameter and a mass 
771 — 118.36 g. The picture of the pendulum in Figure [5^ also exhibits a paper clip and 
the nylon thread, the latter being solidly tied to the brass ball and the paper clip. 

The paper clip is secured to the hook of the chosen counterweight, in turn picked out 
from nine hook masses whose measured values are M = lO.Olg, 20.02g(x2), 50.05g, 
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Figure 4: Photo of the Swinging Atwood Machine (SAM): a pendulum (on the left) and a counterweight 
(on the right), linked together by a nylon thread. 

lOO.lOg, 200.22g(x2), 500. 51g, and lOOO.lOg (Figure [St). The relative difference be- 
tween tfiese values and those engraved in each hook mass is 0.1%; thus, with respect 
to the orders of magnitude of the different masses involved in the experimental device, 
this difference can be neglected. The values considered are therefore presumed to be 
those indicated on the hook mass themselves, namely M — lOg, 20g(x2), 50g, lOOg, 
200g(x2), 500 g, and 1000 g. Henceforth, and for the sake of linguistic simplicity, these 
hook masses will be called "counterweights" , although weight and mass are different no- 
tions, however related. This set enables varying the counterweight mass from 10 g (one 
mass) to 2 100 g (addition of all the masses) with a step of 10 g by hooking several masses 
together. Among these hook masses, one is hung on the nylon thread by means of the 
paper clip, whose measured mass is 0.37 g. It is interesting to note that, by a fortunate 
coincidence, the mass of the paper clip is equal, with a 0.01 g difference, to 0.36 g, i.e. 
the mass of the brass ball minus 118 g. Therefore, the mass of the brass ball can be taken 
as equal to 118 g and the mass of the paper clip can be ignored. Finally, we get, for the 
pendulum and the counterweight, respectively: m = 118 g and lOg < M < 2100 g. 

4-2. About the nylon thread and the pulleys 

The thread (Figure ^) ensures a mechanical coupling between the pendulum and 
the counterweight through the two pulleys. The length of the thread is about one meter 
and its measured mass of 0.10 g is negligible compared to the other masses involved. In 
addition, the nylon thread is assumed inextensible. During experimentation, no thread 
breaking has been reported. 

Pulleys used are shown in photos of Figure [6l They are made up of two parts: 
an internal, immobile one bound to the revolution axis, and a mobile, external one 
liable to rotate around this axis. These two pulley components are uncoupled through a 
ball bearing which, moreover, reduces mechanical energy dissipation by friction. Pulley 
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(a) 



(b) 



Figure 5; (a) Brass ball (pendulum), nylon thread and paper clip (b) Set of hook masses used in the 
experiment (counterweight). Prom left to right: 10 g, 20 g, 20 g, 50 g, 100 g, 200 g, 200 g, 500 g, 1000 g. 

radius is R — 2.5 cm and that of the motionless part is 1 cm. Pulley V, associated to the 
pendulum, has been modified in order to make its groove deeper. Indeed, during the first 
experimentations we observed that the thread could rapidly exit the groove because of 
the pendulum motion. To avoid this, which could by the way be dangerous, two metallic 
plates were added and fixed to the pulley in order to increase by 1 cm the depth of the 
groove (Figures [6^ and \6jp). It is worth noting that the plates are fixed to the immobile 
part of the pulley and are in no way in contact with the mobile one. The motion of the 
latter one is thus not affected by such a modification: hence, from a mechanical point of 
view, the resulting pulley is identical to the original one. 

4-3. Strengthening of the apparatus 

Figure 2] also features two horizontal metallic rods binding the two feet of SAM. 
Their role is to reinforce the machine. Indeed, due to considerable stress involved in the 
pendulum and counterweight motions, the orientation of the two pulleys, as well as the 
arbitrary distance D = 57.5 cm between them, can change; thus, it could be dangerous 
not to strengthen the whole device. The lower metallic rod is solidly fixed to the feet 
while the upper one is solidly fixed to the revolution axes of the pulleys. Consequently, 
we ensured a constant distance between the pulleys whose axes keep a constant direction; 
in virtue of such a construction, SAM is solid and operational. 

5. Motion of the pendulum: experimental results 

5.1. Experimental measure of Ip 

The presence of the inertial momentum Ip of a pulley in ([5]) renders an experimental 
determination thereof necessary. This measure was made using a simple Atwood machine 
where the heavier mass (M = 130 g) fell down from a convenient height h = 1.66 m; the 
lighter mass (the brass ball) being m = 118g. Using equation ([7]), one gets: 



= R 



,2 UM -m)g{At)^ 
2h 



(M + to) 
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(c) (d) 

Figure 6: Photos of the two pulleys, (a) Pulley V, front view: a metallic plate is added which hides 
the internal part and the ball bearing of the original pulley, (b) Pulley V, top view: the groove is made 
deeper by 1 cm with the additional plate; the original pulley is easily recognizable between the two plates, 
(c) Pulley V' , front view: the ball bearing uncouples the immobile internal part and the mobile external 
one (d) Pulley V' , top view: photograph of the axis of revolution and the groove of the pulley. Moreover, 
the two pulleys are distant enough from one another to avoid a pendulum-counterweight collision during 
the motion. 

where At is the faU duration. Through a set of ten measures with a chronometer of 
0.01s of accuracy, the mean fall duration found is (At) — 2.70 ± 0.01s. (Ip) ~ 6.85 x 
10~^kg • m^ ensues. Concerning the value of Ip, errors are due to the measure of both 
At and the positions of M at the initial and final times ~ that is, the determination 
of h. Uncertainties are mainly due to the determination of the final time, which must 
correspond to the falling distance h as precisely as possible; initial and final positions are 
determined with an error of 0.1 cm, which compared to the value of h can be neglected. 
Consequently, uncertainties in position determinations are disregarded and the error 
in the measure of Ip can be reasonably associated to the uncertainty in At (0.4 %). 
Thus, precision on Ip is twice that of At, hence about 1 %; absolute uncertainty is thus 
0.07 X 10~^kg • m^. Therefore, we can write: 

/p = (6.85 ±0.07) X 10"^ kg- m^ or /p = 6.85 x 10"^ kg • m^ ± 1 %. 

5.2. Experimental results 

The motion of the pendulum has been filmed for various /i- values and initial conditions 
(''0 ; ^o)- Then, using the "Synchronie" software and focusing on each film image by 
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image, a pointer enabled us to pick up pendulum positions and record them. Such a 
process is necessarily a source of errors, as it is sometimes difficult to locate the pendulum 
exactly, especially if velocity is high. The errors introduced by such a procedure are not 
simple to estimate. However, trajectories have been correctly recorded, as comparisons 
with numerically-simulated theoretical results will show (see Section . 

5.2.1. Case /xthco = 3 

Since previous studies focused mainly on the particular and theoretical case /xthco = 3, 
this was the first one we experimentally addressed. In fact, the masses available only 
allowed us to approach ^itheo- with m — 118 g and M — 350 g one obtains /ioxp — 
2.966, which is the closest value to /ithco- The sampling time step has been 67ms. 
Motion has been researched for four different initial conditions (ro ; 6'o)=(0.649; 53.5), 
(0.710; 66.5), (0.854; 68.3) and (0.867; 51.1); vq is in meters, 9q in degrees. The motion 
of the pendulum presents the same pattern and characteristics for all these conditions, 
so only the trajectory for the first initial conditions is shown in Figure [7^. All in all, 
359 sampling times have been recorded. The pendulum has a planar revolving trajectory 
around the pulley and presents an asymmetry with respect to the vertical direction. 
Note that the pendulum becomes closer and closer to the pulley as a consequence of 
dissipative phenomena and is bound to end up knocking against it. Phenomena qualifying 
as dissipative are, to our knowledge: the friction between the thread and the pulley, 
the air friction on the pendulum and the counterweight as well as friction inside the 
ball bearing. Evolutions of the length of the pendulum r and angle are displayed in 
Figure [7|d and Figure [71: respectively. The asymmetry of the trajectory and dissipation 
are observable in the evolution of r since this variable exhibits different minimal and 
maximal amplitudes which decrease in function of time t. 

The Fourier analysis (non-displayed) of the data shows that, for the behaviour of 
9, the most relevant harmonic is the constant term, followed by harmonics 7, 6 and 
8. They account, respectively, for 0.401, 0.311, 0.117, and 0.091 of the total variation. 
Of course, the contributions of harmonics 6 and 8 are due to leakage of the 7**^ one. 
The amplitude of the 6'^ harmonic is larger than that of the 8**^, showing that the 
true dominant average frequency is slightly less than 7 times the basic frequency. From 
the number of data and the sampling step time, a basic frequency of 0.042 Hz follows. 
Hence, an average value of the dominant frequency can be estimated to be equal to 
0.294 Hz. However, as is clear by looking at the maxima of the angle, the frequency 
changes with time. The spacing between successive maxima takes the approximate val- 
ues 3.820s, 3.909s, 3.379s, 3.311s,3. 173s, 2.780s, and 2.643s. Hence, the instantaneous 
frequency changes from about 0.262 Hz to about 0.378 Hz. The explanation is simple: 
the dissipation reduces the energy and the length of the pendulum becomes shorter, 
increasing the frequency. 

For the radius, the major contribution comes from the constant term, followed by 
harmonics 13 and 7. They account for 0.810, 0.037, and 0.026 of the total variation. 
Comparing the plot of r as a function of time with that of 0, one can see a doubling 
in the number of maxima. The explanation is simple: largest maxima occur at the left 
part of the plot of the orbit. Then, a minimmn is reached when 9 = — tt (i.e., upwards), 
followed by a maximum to the right, and a new minimum at = —it to reach a larger 
maximum on the left. The spacing between successive larger maxima of r is very close to 
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the one observed for 9. As mentioned, the largest non-constant harmonics are the 13'^ 
and the 7*'', not in a 2-to-l ratio. This is related to the fact that the "best" estimate 
of the main frequency for 6 is slightly less than 7 times the basic frequency. The closest 
integer to the double would be 13 rather than 14. 

The decrease in energy can be calculated as follows. From the experimental data, the 
values of r and 9 can be computed. For that, we have used two independent methods. 
The first one is simply numerical differentiation with a central formula. The second 
one aims at filtering errors in the data as well. A Discrete Fourier Transform has been 
computed and harmonics up to order 32 have been retained. Then, it is possible to check 
that the reconstruction agrees quite well with the initial data and one can compute the 
values of r and 9 using these Fourier expansions. To prevent leakage due to the fact that 
the data at the ends of the interval are quite different (Figs. [7)3 and c), which originates 
the well-known 0{l/n) decrease of the order of magnitude of the n'^ harmonic, different 
procedures have been used, but the results are essentially the same. They also show a 
reasonable agreement using the first and second methods. 
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Figure 7: Case fiexp = 2.966 with initial conditions ro = 0.649 m and 6o = 53.5°. (a) Experimental 
positions (black crosses) and interpolated pendulum trajectory (solid line). The Cartesian coordinates 
of the initial position are y = —0.54 m and z = 0.37 m and those of the final position are y = 0.02 m and 
z = —0.05 m. The pulley is represented by the circle whose centre is at the origin of coordinates (0; 0). 
(b) Experimental positions (black crosses) and interpolated curve (solid line) for the evolution of r. (c) 
Same as (b) but for the angle 6. 
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When f and 9 are available, one can compute pr and pe and subsequently the value 
of the energy. The values for which reaches a maximum (i.e., on the left of Figure [7^) 
are shown in Figure [51 The rate of decrease of the energy is about 0.037 J • . 



s 




Figure 8: 
We have 



5 10 15 20 25 

i(s) 

Values of the energy (crosses) occurring when d reaches a maximum at the left part of Figure[7^. 
also displayed the function Sm = 1-723 — 0.037 1. 



5.2.2. Case /Utheo = 1-5 

The experimental value of /i closest to /ithco = 1.5 is /ioxp = 1.525, given by M = 180 g. 
Proceeding as in the above case allows us to retrieve the experimental trajectory and 
the evolution of the degrees of freedom. In this case, time step is 40 ms. Two initial 
conditions have been considered: (ro ; 6'o) = (0.484 m ; 87.0°) and (0.621m; 87.7°). Since 
they produce the same dynamic behaviour, only the first one is displayed (Figures [5^, 
andint respectively). A slight asymmetric trajectory with respect to the vertical direction 
and an evolution of 6 close to periodic with a period around 1.1 — 1.2 s can be observed. 
For the evolution of r, asymmetry and slight dissipation are also observed. 

A study similar to that of /ioxp = 2.966 is performed. The main contribution to the 
Fourier analysis of 9 comes from harmonic number 7 which accounts for 0.912 of the 
total variation with a frequency of approximately 0.875 Hz, although a better value for 
the average frequency seems to be 0.89 Hz. Harmonics 14 and 21 also play a relevant 
role. As we did for /icxp = 2.966, we can consider the spacing between successive maxima 
of 9 which takes the values 1.125 s, 1.119 s, 1.114 s, 1.110 s, 1.156 s, 1.133 s and 1.063 s, 
showing a decreasing trend with irregularities. 

For r, the largest harmonic is the constant term which accounts for 0.971 of the signal. 
If we skip this term, harmonics 7, 14 and 3 are clearly seen. They contribute to 0.30, 
0.29, and 0.14 of the signal minus the constant part. 
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Figure 9: Same as Figure[7]but for /^exp = 1.525 and initial conditions rg = 0.484m and Bq = 87.0°. In 
(a), the Cartesian coordinates of the initial position are y = —0.485 m and 2 = 0; for the final position, 
one has y = —0.330 m and z = 0.195 m. The origin (0; 0) corresponds to the centre of the pulley (not 
represented). 
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The decrease of the energy as a function of time has been displayed in Figure [TOl this 
time using the values of the energy computed at the minima of 9, on the right part of 
Figure [nti- Now the rate of decrease is about 0.024 J s^^. In this case, using the filtered 
Fourier methods gives better results, because of the large changes in position with a time 
step of 40 ms. 




Figure 10: Values of the energy (crosses) occurring when 9 reaches a minimum at the right part of 
Figure[9^. We have also displayed the function £rn = 0.751 — 0.024 



5.2.3. An unbounded case: /ithco = 1 

In this situation, the experimental value of fji is /Zoxp = 1-017 and three initial con- 
ditions have been considered: rg = 0.120 m, 0.263 m and 0.477 m for Oq = 77.9°. Again, 
trajectories present the same pattern, so only one is shown (Figure [TTkl. They are 
characterized by an increase in r and 9 oscillations around the vertical {9 = 0) with a 
decreasing amplitude of 9 (Figures lllb and c respectively). In Figure [TTb. r appears to 
approach a linear increase in time: r = Vzt where Vz is the velocity along the vertical. 
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Figure 11: Same as Figure[7]but for /^exp = 1-017 and initial conditions rg = 0.120m and 9o = 77.9°. In 
(a), the Cartesian coordinates of the initial position are y = —0.122 m and 2 = 0; for the final position, 
one has y = —0.116 m and z = 1.167 m. The origin (0; 0) corresponds to the centre of the pulley (not 
represented). 
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6. Numerical solution of SAM equations of motion 

6.1. Theoretical trajectories 

Equations of motion Q have been numerically integrated for the same initial condi- 
tions and values of the parameter fi as above in order to compare the theoretical trajecto- 
ries, displayed in Figured^ with the experimental ones. Computation of the theoretical 
trajectories has been performed by using both a Taylor integration method and a vari- 
ety of Runge-Kutta methods of different orders with step-size control to integrate the 
equations of motion. 

From a general point of view, the theoretical trajectories seem quite similar to the 
experimental ones. However, some slight differences can be detected. Firstly, it is obvious 
that dissipative phenomena, though experimentally reduced, play a non-negligible role 
since convergence of the pendulum towards the pulley for the first case (Fig.[7)D), decrease 
of r for the second one (Fig. ^jp) and relatively slow increasing of r for the third one 
(Fig. [Tib) are clearly associated to energy dissipation. Friction will be studied a bit 
further in Subsection [513 and much more in future works. Secondly, it must be said that 
for the case /Xoxp = 2.966 the pendulum touches the nylon thread at each revolution, an 
effect which could be included into the equations of motion through a dissipative term; 
however, it seems quite difficult to introduce such an effect in a realistic manner. 

An important point concerns the comparison between these trajectories and those of 
(sot - For the first case for instance, Tufillaro's trajectories are symmetrical with respect 
to the vertical direction, as opposed to the above ones (Figures [7^ and [T2k). Clearly, 
this asymmetry is due to the influence of the pulleys. 

6.2. Influence of the pulleys on the motion 

Pulleys can influence the motion of SAM through their dimension (since radius R ^ 0) 
and their rotation (since Ip ^ 0). 

6.2.1. Influence oj the radius: /icxp = 2.966 

Figure [T3| sketches the trajectory for /Zexp = 2.966 when Ip — Q for different values of 
increasing i?; i.e. the pendulum pulley V has a non-negligible radius but pulleys are not 
allowed to rotate. The first figure corresponds to the symmetrical Tufillaro trajectory 
[R = 0, /p = 0). Ostensibly, the larger R is, the more significant the asymmetry becomes; 
for R = 5 cm, the pendulum hits the pulley before completing one revolution. 

6.2.2. Influence oJ the inertial momentum: /ioxp = 2.966 

Figure [T4| sketches the trajectory for /ioxp = 2.966 for increasing values of Ip with 
a value of pulleys radius taken to the real value R — 2.5 cm. The trajectory of the 
pendulum is visibly modified: it describes more irregular trajectories and fills more space 
as increases. 
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Figure 12: Numerical trajectories obtained by solving equations of motion ((5J. (a) ^icxp = 2.966 and 
for the same initial conditions as Figure [7ji. (b) fiexp = 1.525 and for the same initial conditions as 
Figure [5^. (c) /.tcxp = 1.017 and for the same initial conditions as Figure [TTJi. 
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Figure 13: Theoretical trajectories of the pendulum with fixed pulleys of different radii R for /Xexp = 2.966 
and initial condition (ro ; eo)=(0.649 m ; 53.5°). (a) R = 0, {h) R = 2.5cm, (c) R = 4cm and (d) 
R = 5 cm. In the latter subplot, the pulley has been represented by the circle whose centre is at the 
origin of coordinates in order to see the collision between the pendulum and the pulley. 
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Figure 14: Theoretical trajectories of the pendulum with pulleys of radius R = 2.5 cm and different 
inertial momenta Ip for /iexp = 2.966 and initial condition (ro ; So)=(0.649m; 55.7°). (a) Ip = 0; this 
subplot which is the same as Figure 113b has been repeated for the sake of clarity, (b) Ip = 6.85 X 
10-'' kg ■ m^, (c) Ip = 13.70 X lO"** kg ■ m^ and (d) Ip = 27.40 X 10"'^ kg ■ m^. 
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Figure 15: Theoretical trajectories of the pendulum with fixed pulleys of different radii R for /Xexp = 1-525 
and initial condition (ro ; eo)=(0.484 m ; 87.0°). (a) _R = 0, (b) _R = 2.5cm, (c) R = 5cm and (d) 
R = 10 cm. 
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Figure 16: Theoretical trajectories of the pendulum with pulleys of radius R = 2.5 cm and different 
inertial momenta Ip for /icxp = 1.525 and initial condition (ro ; 6o)=(0.484m; 87.0°). (a) Ip = 0; this 
subplot which is the same as Figure 115b has been repeated for the sake of clarity, (b) Ip = 6.85 X 
lO-** kg ■ m^, (c) Ip = 13.70 X lO-** kg ■ m^ and (d) Ip = 54.80 X 10"'^ kg ■ m^. 
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6.2.3. Influence of the pulleys: /^cxp = 1-525 

In this case, results of variations in R with Ip — Q (Figure llSp and Ip variations for 
R — 2.5 cm f Figure [TB]) are similar in that trajectories are modified more visibly as R and 
Ip increase. However, the influence of Ip and R on SAM motion depends on the value 
of /X considered. If R increases, the pulleys being fixed, the brass ball ends up hitting 
the pulley but asymmetry does not becomes more and more important. If Ip increases, 
pulley dimensions being fixed, the pendulum evolves in a much more limited space. 

6.3. Poincare maps and rotation number 

To have a global view of the dynamics of SAM, we have computed Poincare maps 
Vm on suitable Poincare sections. We note that the Hamiltonian ^ is not 27r-periodic 
because of the linear term in 9. Be that as it may, our Poincare section S defining Vm 
is given by the coiling of 9 through multiples of 2tt, with ^ > 0. But, contrarily to SAM 
without pulleys ([13]), one has to distinguish between cuts through different mifltiples of 
2iT. Plus, one can not superimpose the different sheets. 

Some types of "escape" can occur. The main source thereof is r going to zero. Other 
relevant sources of escape are r increasing too much or \9\ becoming too large. All orbits 
leading to some of these escapes are deleted. 

To compare with the experiments, we present some examples for M = 350 g and M = 
180 g. As levels of energy, one has taken the values corresponding to the experiments de- 
scribed in[5Xl]and[5221 that is, (ro; 6'o)=(0.649 m ; 53.5°) and (ro ; 6'o)=(0.484m ; 87.0°), 
respectively, with zero initial velocity. Given r, pr in E and 9^0, the value of pg is 
recovered from the energy level. Figure [T7] shows some results. All these massive com- 
putations use Taylor integration methods, in order to ensure a very good conservation 



The top left plot corresponds to M = 350 g, leaving S with pr > 0. The points on 
S with Pr < correspond to cuts through 9 = — 27r, in agreement with the description 
of motion in 15.2.11 To produce the Poincare map, we first computed the periodic orbit 
as a fixed point in E n {pr > 0}. The approximate values (r* ; p*) of the fixed point are 
(0.332814 m ; 0.554330 kg • m • s^^). This periodic orbit can also be obtained by starting 
the motion from rest at (rg ; 9o) w (0.61316m; 64.032°), not too far from the values 
used in the experiment. The curve drawn with large dots around the periodic orbit 
shows the iterates of Vm corresponding to the data used in the experiment. It is clear 
that the theoretical, non-dissipative, motion seems to be in a 2D torus, but numerical 
computations can never exclude the possibility of having a periodic orbit with very long 
period or a tiny chaotic zone. The intersection of this torus with 9 — ~2n is also shown 
in the region pr < 0. The asymmetry is ostensible in the plot, and is due to the effect of 
the pulleys. To produce the full plot, we have taken initial conditions on E with r = r* 
and different values of Pr starting at p*. From some value of Pr onwards, iterates escape. 
We do not exclude the presence of tiny islands outside the last invariant curve shown. 

For the sake of completeness, we show in the top right plot the Poincare iterates 
leaving S in the region {pr < 0}. Points appearing in {pr > 0} are on the sheet {9 — 27t}. 
The fixed point in {pr < 0} is approximately (0.383367m; —0.620020kg- m - s^^)- 



of the energy (see, e-g.. 




)■ 
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Figure 17: Iterates of Poincare maps for M = 350 g (top row) and M = 180 g (bottom plot) on the 
levels of energy of the physical experiments described in the text. The variables (r ; pr), as usual in m 
and kg ■ m ■ s~^, respectively, are shown for 9 = (mod 2tt) and 8 > 0. On the top left plot, the points 
leave the section with pr > and = 0, and those with pr < are in 8 = —2tt. On the right top row, 
top initial points are in 6 = with pr < 0, and those showing in pr > have 9 = 2it. On the bottom 
plot, all intersections occur at = 0. The thick curves are the intersections of the theoretical 2D tori 
corresponding to the physical experiments with the section S. 



The bottom plot corresponds to M — 180 g. In that case, only intersections having 
6* = are found. The fixed point is (r* ; p;) w (0.356386 ni ; 0.008848 kg • ni • s^^), which 
can also be obtained leaving from (rp ; ^o) = (0.53313 m; 78.596°), again not too far 
from the values used in the experiment. As before, the curve drawn with large points 
around the periodic orbit would be the one obtained for the physical experiment without 
dissipation and, as expected, denotes motion on a 2D torus. 

A useful tool to understand the dynamics of Area-Preserving Maps and, in particular, 
Poincare maps such as the ones displayed, is the rotation number p for the map restricted 
to invariant curves. Despite the fact that the rotation number still exists for periodic 
orbits of Vm and for the eventual islands around them (thereupon being rational), it 
is not defined, in general, for orbits with chaotic dynamics. The method used for the 
computation is topological and based on the order of the arguments of the iterates with 
respect to the central fixed point (r* ; p*) of Vm- The procedure computes two estimates 
Pinf and psup such that pinf < p < Psup- If for some orbit one has pinf > Psup this proves 
p not defined. If the number of Vm iterates is iV, the typical errors, when p exists, are 
O (A^^^) for constant type rotation numbers. See the Appendix in '23'| for details and a 
complete analysis of the error depending on the Diophantine properties of p. 
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In Figure [TSl we show results corresponding to the Poincare maps displayed in Figure 
1171 top and bottom. The computations are done starting at initial points of the form 
(r* ; Pr) with ^ P* ~ J A, j — 1,2,... with a small step A. As successive iterates fall 
in E for values of alternating between and — 27r, the map has been used instead 
of Vm- On the left plot, we display p as a function of Pr for M = 350 g. This is the curve 
which has a large dot near the upper right corner. The point corresponds to p^ = p* and 
the limit rotation number. We see a decreasing rotation number when pr decreases up to 
a value pr ~ 0.143982kg • m • s~^. Beyond that point, the Vm iterates escape. The other 
curve, also shown here for comparison, corresponds to the pulley-less case with ^ — 3 
and will be described later. 

In fact, what seems a nice curve for M = 350 g should have, generically, a "devil's 
staircase" structure. That is, there are infinitely many intervals in which p € Q and it is 
constant. They correspond to islands around elliptic fixed points. Some of these rational 
values, such as 1/6, 1/7, 1/9, 3/20, . . . (or resonances) have been detected. But they are 
very narrow. As an example, the inset in the left plot of Figure [TSl shows the behaviour 
of p in an interval of pr whose width is lower than 10^'^ kg • m • s^^ and p around 1/9, 
which illustrates a typical pattern when crossing a resonance through an island. 
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Figure 18: Plots of rotation number for different cases. Left plot: the rotation number p as a function of 
the initial value of pr with r = r* for M = 350 g and also for the pulley-less case with = 3. A zoom-in 
around a resonance is shown in the inset. Right plot: p as a function of pr for M = 180 g. See the text 
for additional details. 



On the right hand side of Figure [T51 we show the results for M = 180 g. As before, 
we use initial points (r* ; pr) with pr going away from the fixed point of Vm- V^^i 
also been used instead of Vm because the latter is close to a symmetry with respect to 
(r* ; p*), hence the displayed rotation number is small. It increases monotonously as we 
move away from the fixed point, also marked as a large dot. Now, however, the intervals 
with p G Q constant are extremely narrow. We observed a few resonances, such as 2/45 
and 5/116, checking that the width of Pr-intervals is below 10~^ kg • m • s'-^. 

We want to mention now that the pattern of p as a function of Pr is a clear indication of 
non-integrability. Indeed, let us first look at the left plot in Figure[T51 We have also shown 
the results for SAM without pulleys and n — 3, the integrable case. The level of energy is 
the one corresponding to (rg ; 9o) = (0.25 m ; 0.0°), similar to the kind of reference orbits 
used in [lj|. A periodic orbit is found near (r* ; p*) = (0.125 m, 0.369587kg • m • s~^). 
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Using r = r* and Pr between p* and zero, we check that p decreases to zero. Under V^^ 
there is a "separatrix" bounded by r = 0, = and a curve of the form r = 0.25 — ap^ 
for a suitable a. The system being integrable, the map V^^ has a first integral / and the 
iterates are on level curves of /. When approaching the separatrix, the dynamics slows 
down near r = and pr , hence p is very small. 

But in the pulley case, the rotation number ceases to exist at a value which is unrelated 
to any separatrix. This is against the typical behaviour of integrable maps. 

One should expect chaotic dynamic regions in a resonance zone in SAM with pulleys. 
Beyond the islands around periodic elliptic points, there are homoclinic tangles associated 
to the hyperbolic zones. Attempts to visualize them lead to the conclusion that the size 
of the "chaotic zones" is, at most, of the order of magnitude of the round-off errors with 
standard double precision computation. Hence, the escape remains a main evidence of 
non-integrability. 

6.4- Remarks on dissipation 

As we have already noted, experiments show non-negligible dissipative phenomena 
which decrease mechanical energy during the motion. Hence, equations of motion such as 
([5]) which do not contain any dissipation term yield by no means a complete description 
of SAM dynamics. 

In order to show that the observed convergence of the pendulum towards the pulley 
(Figure [7^) is due to dissipation and not to a too short integration time interval, the 
corresponding theoretical trajectory (Figure \T^ ) has been integrated for much larger 
intervals. Trajectories are not shown but we can say that the limited region of space 
occupied by the pendulum in Figure 112b is progressively filled as time is running. The 
corresponding Tufillaro trajectory integrated for larger time is in addition symmetrical. 

Four main different sources of dissipation can be noted: air friction on the pendulum 
and the counterweight, thread friction on the pulley grooves due to an imperfect absence 
of slippage condition, and dissipation due to the ball bearing of the pulleys. This last 
source is implicitly included in the equation of motion through the measure of Ip since 
it has been determined using a pulley (see Sectional). 

An estimation of air friction can be performed by comparing the weight of the pen- 
dulum, mg, to the drag force exerted by air on the pendulum, whose expression is: 

p _ ^ PaV^ o 
i'd — L-rf 2 J 

where Cd is the drag factor, pa — 1.29 kg • m^'^ the air density at ambient temperature 
T = 298 K, V the velocity of the ball, and S the effective surface of the pendulum. For 
a spherical ball, Cd ~ 0.4. Since v ~ Im-s""'^, S — TrD^/A with Df, = 30mni, the 
drag-weight ratio is of the order of 2 x 10^"*. One gets the same order of magnitude for 
the cylindrical counterweight for which ~ 1. Dissipation will be the detailed topic of 
a further paper. 
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7. Non-integrability of SAM with pulleys 

In this Section, a rigorous and original analytical proof of non-integrability of SAM is 
performed in order to complete the above rotation number and Poincare section analysis. 
Needless to say, this non-integrability result is fundamentally different, both in approach 
and scope, from the numerical and graphical inference. 

Since SAM may be expressed in terms of symplectic formalism, in order to detect or 
predict chaotic behaviour, it is pertinent to recall some concepts related to the integra- 
bility of Hamiltonian systems in the sense of Liouvillc- Arnold. 

7.1. Algebraic background for studying integrability 
7.1.1. Linear and Hamiltonian integrability 

Differential Galois Theory.. See [l5l |. (s^ and [lol | for more information. Given a linear 
differential system, with coefficients in a differential field {K, d) whose field of constants 
C is algebraically closed (e.g. [C{t) , d/dt]), 

dy = A{t)y, (8) 

an algebraic group G exists, called the differential Galois group of ([8]), acting over the 
C- vector space {ipi, . . . , of solutions of ([8]) as a linear transformation group over C. 
Furthermore, G contains the monodromy group of ([5]). The Galoisian formalism proves 
useful here due to the following: ^ is integrable (i.e. its general solution can be written 
as a finite sequence of quadratures, exponentials, and algebraic functions) if, and only if, 
the identity component G'^ of the differential Galois group G of ^ is solvable. 

Everything said in the previous paragraph may be obtained, mutatis mutandis, for 
linear homogeneous differential equations 

d" d""i d 
an {t) -^y + a„-i {t) -^j^y H '^ai{t)—y + aa (t) y = 0. 

In Subsection 17.1.31 we will denote Galois groups in this Gal (L) accordingly, L S 
C{t) [d/dt] being the corresponding differential operator. See [l2| for more details. 

Liouville- Arnold integrability.. On the other hand, we call a Hamiltonian system q — 
dH/dp, p = ^dH/dq, whether or not linear, meromorphically integrable (in the sense of 
Liouville- Arnold) , if it has as many independent meromorphic first integrals in pairwise 
involution as degrees of freedom. Same applies in the above definition, mutatis mutandis, 
substituting algebraic, rational or any other function class in for meromorphic. For 
the sake of simplicity, conjugate canonical variables will be henceforth assembled in a 
single vector z — {q,p) and the Hamiltonian system will be written in compact form 
z = Xh (z). Everything is considered in the complex analytical setting from this point 



7.1.2. Morales- Ramis-Ziglin Theory 

For each integral curve F = : t G /| of a given autonomous dynamic system in 

dimension m 

z^Xiz), (9) 
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the variational equations of order k for ([9]) along F, VEp, are satisfied by 2^ := d^cfi/dz^ 
- see, e.g., [l7| for explicit expressions of VEp for general k in terms of vectors k — 
(fci, . . . ,km) G such that k ~ fci + . . . + km- Note that d^c^jdz^ is an abridged 
notation for the {k + l)-dimensional matrix of all partial derivatives of 0. In other 
words, it's a vector for fc = 0, a matrix for fc = 1, etc. 

In particular, for fc = 1 and denoting the matrix of the first-order variational equations 
simply by S, we obtain 

S = X' (0) S. (VEr) 

We thus have: 

• a (generally nonlinear) system ^ and 

• a linear system dVErD linked to ([9]). 

The hallmark theorem in this approach connects the two notions of solvability listed in 
17.1.11 namely as applied to a Hamiltonian Xh and its variational equations VEr, along 
an integral curve F of Xh- The whole theory is actually the ad hoc implementation 
of the following heuristic principle: if a Hamiltonian is integrable, then its variational 
equations must also be integrable. 

We assume F, a Riemann surface, may be locally parameterized in a disc 7 C C of 
the complex plane; we may now complete F to a new Riemann surface F, as detailed in 
[igL §2.1] (see also [l5l, §2.3]), by adding equihbrium points, singularities of the vector 
field and possibly t — 00. 

Theorem 7.1 (J Morales-Ruiz & J-P Ramis, 2001). Let H be an n- degree- of- freedom 
Hamiltonian having n independent first integrals in pairwise involution, defined on a 
neighbourhood of an integral curve F. Then, the identity component Gal (VEp)*' is an 
abelian group (i.e. Gal (VE)p) is virtually abelianj. 

See [igI . Corollary 8] or [l^ Theorem 4.1] for a precise statement and a proof. 

7.1.3. Differential operators. A primer in the Boucher-Weil Theorem 
Linear differential equations.. See [26l §2] for more details. Let 

/d\" fdV'-' 
-^ = "" 37 +""-1 37 ^ ^"0; a„,...,ao e L(i) 



dtj ""^ Vdi, 

be a differential operator with coefficients in the field of formal Laurent series. If L (y) — Q 
has a solution of the form y — t^ J2k>o '^'^^^ ^^"-^ '^o 7^ 0; the formal substitution of y 
into the differential equation yields the vanishing of all powers of i, the smallest one 
among them - we call the equation P (p) ~ derived from the latter vanishing the 
indicial equation (at 0), the roots of which are usually called exponents of L(ij) (at 
0). In particular we can also define the indicial equation at infinity by means of the 
transformation / = 1/t and expansion around / = 0. 

It is a known fact ([2^ Lemma 2.1]) that the degree of P {p) = is at most n. A 
singular point c of L is called regular singular if deg P (p) = n. A linear differential 
equation L = with only regular singular points (including 00) is called Fuchsian. 
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We call L reducible if it factors in a product of operators of smaller positive order. 
An operator L admits a first order factor d/dt — /, / G C (t) if and only if L{y) =0 
admits a solution y such that y ~ fy; in particular: 

Lemma 7.2. If L is of order 2: L is reducible if, and only if it has an exponential 
solution, i.e. a solution whose logarithmic derivative is rational. □ 

Lemma 7.3 ( [2^ . §3.1.2]). If L is Fuchsian, every exponential solution must be of the 
form y = Yli^iit — tiY^ P (t) , where P G C [X] andti,...,ts are finite singularities of L 
with exponents ei, . . . ,es, whether or not integers. 

Normal variational equations. Let 

Id„ 



' -Id„ 

be the canonical symplectic matrix. Given a Hamiltonian system i — Xh {z) = JWH (z) 
expressed in Darboux canonical coordinates z = {q, p) , system dVErD reads S = JH" {z) E 
along T = {z (t)}. 



Gauge transformations.. A P [A] .= P ^ \^AP — Pj, P being a given symplectic 

matrix, may be used to reduce dVErD by selectively vanishing a number of entries in 
JH" {z) (see 0, §5.2], [3, §4-1], (3|). A typical first choice consists of symplectic matrices 
of the form 

P{t)=( ^~Z C2 C3 C4 V (10) 



in order to induce a row and a column of zeroes in the variational matrix. When such is 
the case, ( |VErD acquires a "decoupled" appearance and a system of order n ~ 2 may be 
extracted therefrom. Such a system is usually called a normal variational system along 
r. Let us denote it by NVEp. 

As is always the case with all differential systems, a cyclic vector (Q, [1^) may be 
used to obtain a linear differential equation L (y) ^ of order n — 2 equivalent to NVEp. 

The result central to this paragraph, and a particularly useful consequence of Theorem 
17.11 is the following (see also 0, Proposition 9 & Theorem 8 (§5.3)]): 



Theorem 7.4 (D Boucher & J- A Weil, [Sj, Criterion 1]). Let Xh be a Hamiltonian 
system and L its normal variational operator along a given particular solution. If L is 
irreducible and displays logarithms in a formal solution, then Gal (L)" is not abelian, i.e. 
Xh is not integrable. □ 

Typically, the most difficult part in trying to apply Theorem 17.41 is to check the 
irreducibility condition. 

7.2. Statement of the main result 

We recall Mt = M + m + 2Ip/R^ and 

' pI , {P2 + Rpi ' 



+ gqi{M — TO cos 92) — gR{Mq2 — m sin 92), 



Mt mqf 

where qi = r, q2 = 9, pi ~ pr, and p2 ~ pe- The main result in this Section is the 
following: 
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Theorem 7.5 (Non-integrability of SAM with massive pulleys). For every phys- 
ically consistent value of the parameters, regardless of Ip and R, X-h is meromorphically 
non-integrable. 

This is a complement to what has already been proved for SAM without pulleys, i.e. 
the limit case Ip = 0, R = and Mt = M + m: 



■N ^ (pi _^ pI 



gqi {M - m cos 92) ; 



11) 



2\Mt ' mql 
in that case, the following held: 

Theorem 7.6 {SAM without massive pulleys). 

1. (d, Theorem 1], Equation (16)], t^; see also Remark mT2l) If M > m and 

M p(p+l) 

M = — 7^ Aip / , , T 

m p[p + 1) — 4 

for every p €z Z,p > 2, then Hamiltonian X-j-c^ is non-integrable. In particular, if 
fi S (3/2, 3) U (3, 00), it is non-integrable. 

2. ([2^, Equation (16)]) For p = 2, fi — 112 = ^, X-h^ is integrable and has the 
following first integral: 



I = qi42 qic 



91 92 



2 2 2 2, PmC^ - 2P2S 

gq^sc = gq^^c s + p2- 



Avn?qi 



where c = cos ((72 /2) , s = sin ((72 /2) . 



3. ([ij, Theorem 4]) The degenerate cases fip, p > 2 referred to in item 1 are non- 
integrable. □ 

The last case is significantly more difficult; it relies on the higher order variational 
equations 17 1 and uses techniques introduced in (isj . 

7.3. Proof of Theorem\T^ 

We have two particular solutions for SAM with massive pulleys: 



g{M ~m) ( t{l - t) 
Mt 

g{M + m) f t{l - t) 



0, 1 - 2t, R(2t - 1; 



2-K 



Mt ' g{M + m) 



l-2t, R{2t - 1; 



System ( |VEr[ ) around F = {zT;^ii{t)} takes the form ^ ~ with 



/ 





l/Mt + R^ai Rai 


\ 








Rai ai 















V 


02 





/ 



32 



where 

AMf , mg'^{M + m)t{t-l) 

Oi — tt: TTTTT, and ao = . 

mg^{M + m)H^{t-l)^ 2Mt 

Notably, dVErD decouples without the need for an additional gauge transformation such 
as ([TU]) as is usual and customary (0, 15|), and as would be the case if F = {2;^_o(i)}- 
See Remark [73] 1 below. System NVEr takes the form $ = _B$ with 

B = 

A necessary condition of integrability is the virtual abelianity (see Theorem I7.ip of 
Gal (NVEr), NVEr, by means of a cyclic vector c and the subsequent gauge transfor- 
mation given by <i> = Q"^*, where Q = { c c + B'^c )^ ([7, §B.4]), will take the 
form 

/ 1 




dt 



2Mt 2{2t- 1) 



\ (M + m)(<-l)t {t-l)t 

in presence of a constant cyclic vector c — (a, 0). The above system is obviously equiva- 
lent to the following hypergeometric ( 1, §15.5]), hence Fuchsian differential operator: 



dV 2(2t- 1) d 2Mt 



dt) t{t-l) dt (M + m)t{t-l)' 

Let us now check the virtual non-commutativity of Gal (L). L has three singularities: 
0, 1, oo. At t = or i = 1, we have local exponents —1 and 0. The formal solution at 
t = has logarithms except in two cases: m = 0, Mt = M + m, both discarded in our 
case since they would correspond to no small mass and no pulley, respectively. Indeed, 
a particular solution is 

~ P [9(M + m)-8Mt]i/^ 3 , [9(M + rn) - 8Mt]i/^ 

yUtj-2^i(^2 2(M + m)i/2 '2+ 2(M + m)i/2 

where 

{a)k{h)k t'' 



2Fi{a,b; c; t) = 



k=0 



(c)fe fc! 



is the Gauss hypergeometric function ([!), §15.1]), {a)k — a{a+l) ■ ■ ■ (a + fc— 1) being the 
Pochhammer symbol. This first solution has a local expansion around < = of the form 

Mt Mt \2 (M + m) + Mt] , ^ , ,x 

yi - 1 + t + ^ + O (t^) . 

M + m 3 (M + mf ^ ' 

We may then obtain an expansion for a second formal solution 



1-^+ 3+^^^-[^^--^(^,+-nt+o(t^)+2^+— 

t M + m \ {M + mf J ^ ^ M + m 
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Keeping the Boucher- Weil Theorem 17.41 in mind, and in presence of the logarithm in y2j 
there is obstruction to integrability if L is irreducible. 

Let us assume it is reducible. In virtue of Lemma 17.21 L admits an exponential 
solution. We recall that L is Fuchsian. The expansion of an exponential solution around 
i = does not contain logarithms, although we have shown there are formal solutions with 
logarithms around the singularities - as has been seen explicitly for t = Q. Thus, those 
without logarithms correspond to the maximal exponents ([3J|), hence the admissible 
exponents at the finite singularities are all 0. Hence, in virtue of Lemma 17. 3[ the only 
possible form for an exponential solution y would be that of a polynomial solution; let 
iV be its degree. Expanding y in increasing powers of t~^, —N is the exponent of the 
leading term, hence an exponent at infinity: Y^^=Q^ki^ — ^k=Q^N~k{^/t)''- 

Now, the exponents at infinity are the two roots of the indicial equation — 3p + 
2Mt/{M + m). Since -N is such a root, this means Mt = -N{N + 3){M + m)/2 with 
N positive. Therefore Mt would be negative, which is physically irrelevant. 

Hence follows that the Swinging Atwood Machine system with massive pulleys is 
always non-integrable with meromorphic first integrals. □ 

Remarks 7.1. 

1. Intriguingly, the solution used for our proof was ZTr,R, which, at least for the pulley- 
less case i? ^ 0, Mt ^ M + m, and although mathematically plausible, has no 
actual physical significance. The solution which would be physically possible for all 
values of R, zq^h, posed further problems with regards to system dVErD and was 
finally discarded in our proof. It is worth mentioning, however, that in the case 



without pulleys, 1J| used precisely the corresponding form of the latter solution 

^0,0- 

The same proof given for Theorem 17.51 may be obtained, analogously, for SAM 
without pulleys Hw, > 1, and the proof yields precisely item 1 in Theorem 17.61 
Indeed, by using the "classical" solution 2:9,0 (which corresponds to the original 
Atwood machine) and an adequate gauge transformation, we obtain the normal 
variational equation: 

2Mt , , (2f - 1) d , , , „ 

This is a Gauss hypergeometric equation. Local solutions are 



and 

{M ~ m) [2tMt + {M - m)] 
2Mt [Mt + (M - m)] 



Int. 



The degenerate cases Mt = and Mt + M — m — Q are of course not physically 
acceptable. The exponents at zero are and 2, the exponents at t = 1 are also 
and 2, and the exponents at infinity are the roots of the polynomial Poo{X) where 
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P^{X) = {M -m)X'^ + 3{M ~m)X -2 Mt. Same as in the proof of Theorem 
17.51 reducibihty would occur only for a polynomial solution, i.e. in presence of an 
integer N such that Poo{N) = 0, implying 

Mt = i (iV + 4) (AT + 1) (M - to) . 

Setting Mt = AI + m (that is, the pulley- less case) we would have an equation for 
H whose solution would be 

_ + 5N + 6 _ {N + 3){N + 2) 
^ ^ + 5N + 2 ^ (A^ + 3)(7V + 2) -4' 

obviously equivalent to the condition in item 1 in Theorem 17.61 for N = p — 2. 
For 11=1, however, the closest thing to such a proof as that sketched in item 2 
is discarding the existence of first integrals with a meromorphic growth at infinity, 
e.g. rational ones. This is due to the fact that the normal variational equation 
around particular solution z (t) = {—gt/2 ; ; —gm ; 0) is a Hamburger equation, 



1 d_ 

t"'"' t dt"'-"' ' dt2 

i.e. an equation with exactly two singularities: a regular one at zero and an irregular 



-7 2/(0-7-2/(0 + 3722/^=0' 



one at infinity [ll|, §17.6]; the solutions are almost Bessel functions: the general 
solution is: 

y{t) = Cit h (2\/2l) + C2t K2 {2V2t^ , Ci, C2 G C, 

where, n ^ given, /„ (z) and Kn {z) are the modified Bessel functions of the first 
and second kind, respectively (0, p. 416], p. 376], [H, p. 185]): 

^ ± f exp[(./2)(t + lA)] ^ . /_„ (.) - /„ (.) 



27ri 7si ' ^ ' 2 sin(n7r) 

both having a branch cut discontinuity in the complex z plane running from —00 to 
(although /„ is regular at 0, whereas Kn has a logarithmic divergence at 0). This 
assures the presence of a non-trivial Stokes multiplier at infinity for the variational 
equation when /i = 1. This implies the following important conclusion: the system 
is not integrable with first integrals which are rational functions ofr, 6, cos 9, sin9, 
Pr, P0, where 9 belongs to a neighbourhood ofO. 

Indeed, any first integral is a function of {q,p), hence must be defined on the 
phase manifold M := C^. We can partially compactify M into M := x C^, 
where stands for the compactification, by means of the inverse stereographic 
projection , of the domain for (qi ; pi) = (r ; pr), whereas the second factor 
C'^ is the {9 ; pe)-plane. Using L- in order to compactify the particular solution 



■ c 



X 



r = {z {t)}, we obtain a Riemann sphere T in P^, whose immersion in 
is contained in {6 = 0}; therefore cos 9, sin9 are holomorphic, hence meromorphic 
on a neighbourhood of F in M. Theorem 17.11 implies the absence of a complete 
set of first integrals which are meromorphic on a neighbourhood of F in the partial 
compactification M. Since any rational function of r,9,Pr,pe,cos9,sm9 must be 
meromorphic in a neighbourhood of F in M, we obtain the last claim in the previous 
paragraph. 
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8. Conclusion and perspectives 

In this paper, experimental and theoretical results concerning the Swinging Atwood 
Machine are presented. Equations of motion with two pulleys are found generalizing 
those of previous studies. An experimental device for SAM is constructed and described 
in detail. Experiments are conducted and the trajectories retrieved from computer video 
analysis closely match those arising from numerical solution of the equations of motion. 
Such comparisons seem to show that the motion is influenced by the non-negligible 
dimension and the rotation of the pulleys; in particular, a non-zero pulley radius leads to 
asymmetric pendulum trajectories. We conclude that pulleys cannot be ignored in SAM 
dynamics. 

Finally, after giving some numerical evidence of the lack of integrability of SAM in 
the absence of dissipation, a complete proof of this fact is given using differential Galois 
theory and the necessary conditions following from Morales-Ramis theory. 

Several perspectives of this work can be considered. First of all, other experiments are 
currently conducted with /ioxp = 2.03 and qualitative preliminary results (not shown in 
this paper) indicate that the dynamics of SAM seems to be irregular - chaotic behaviour 
is expected. Detailed research will be performed in a future work. Secondly, the influence 
of dissipative phenomena on SAM dynamics should be studied; a possible procedure in 
such direction is increasing air friction by coating the pendulum with different materials, 
judiciously chosen so as to induce changes in drag force. Another method could consist 
in placing the counterweight in media more viscous than air, such as water or glycerin, 
and forcing it to evolve therein. Perhaps an ultra fast camera of about 1000 images per 
second could be necessary to pick up much more points and achieve a better resolution 
in pendulum trajectories. 

Concerning the integrability of SAM without pulleys, that is Hamiltonian Hw in (fTT|) 
obtained from ^ by skipping the contributions of R and Jp, the results summarized in 
Theorem 17.61 along with the result shown in the previous section for fi = 1, close the 
problem: the pulley-less case is non-integrable for all values of /x ^ 3. 
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